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ABSTRACT 

We use a ID thermochemical and photochemical kinetics model to predict the 
disequilibrium stratospheric chemistries of warm and hot Jupiters (800 <T< 1200 
K). Thermal chemistry and vertical mixing are generally more important than 
photochemistry. At 1200 K, methane is oxidized to CO and C0 2 by OH radicals 
from thermal decomposition of water. At T < 1000 K, methane is reactive but 
stable enough to reach the stratosphere, while water is stable enough that OH 
levels are suppressed by reaction with H 2 . These trends raise the effective C/O 
ratio in the reacting gases above unity. Reduced products such as ethylene, acety- 
lene, and hydrogen cyanide become abundant; further polymerization should lead 
to formation of PAHs (Poly-Aromatic Hydrocarbons) and soots. Parallel shifts 
are seen in the sulfur chemistry, with CS and CS2 displacing S2 and HS as the 
interesting disequilibrium products. Although lower temperature is a leading 
factor favoring hydrocarbons, higher rates of vertical mixing, lower metallicities, 
and lower incident UV radiation also favor organic synthesis. Acetylene (the 
first step toward PAHs) formation is especially favored by high eddy diffusion 
coefficients K zz > 10 10 cm 2 /s. In most cases planetary compositions inferred 
from transit observations will differ markedly from those inferred from reflected 
or emitted light from the same planet. The peculiar properties of HD 189733b 
compared to other hot Jupiters — a broadband blue haze, little sign of Na or K, 
and hints of low metallicity — can be explained by an organic haze if the planet 
is cool enough. Whether this interpretation applies to HD 189733b itself, many 
organic-rich warm Jupiters are sure to be discovered in the near future. 

Subject headings: planetary systems — stars: individual(HD 189733) 
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1. Introduction 

The most easily studied transiting hot Jupiters are HD 209458b and HD 189733b 
(Fortney et al 2009). HD 209458b is typical of known hot Jupiters (Burrows et al. 2008). 
Transit observations show a visible spectrum dominated by the 590 nm sodium resonance 
lines (Charbonneau et al 2002), the visibility of which implies a clean, haze-free atmosphere 
(Sing et al 2008). HD 189733b is atypical, although its details are controversial (Fortney 
et al 2009). In visible light HD 189733b presents a featureless blue haze with little sign of 
the expected Na and K lines (Pont et al 2008). Lecavalier des Etangs et al (2008) suggest 
that these lines are obscured by enstatite grains mixed up to microbar levels. Swain et al 
(2008, 2009) fit transit data between 1.5 and 2.5 /im with CO, C0 2 , H 2 0, and CH 4 , all 
with distinctly subsolar abundances corresponding to a metallicity of [M/H] ~ —0.5. They 
suggested that adding C 2 H 2 , C 2 H 4 , or possibly NH 3 , would improve their model's fit. On 
the other hand Fortney et al (2009) fit Desert et al's (2009) dayside spectrum at 4.5 /im 
with supersolar CO and C0 2 corresponding to 0.5 < [M/H] < 1.5. 

An organic haze offers another way to explain the appearance of HD 189733b. Here 
we use a photochemical model to address the threshold for organic haze formation in hot 
or warm Jupiters. In Zahnle et al (2009) we noted that temperatures below 1200 K change 
the kinds of chemistry that take place. In hot atmospheres methane is efficiently oxidized 
to CO. Cooler atmospheres are functionally more reduced because H 2 is more stable: 
organic molecules are favored over CO, PAHs (polycyclic aromatic hydrocarbons) and soots 
(disorganized agglomerations of bigger PAHs) might form and precipitate. The PAHs and 
soots would be the source of the haze. 
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2. Previous Models 

The possibility that organic hazes might be important in irradiated brown dwarfs was 
first raised by Griffith et al (1998). However, equilibrium chemistry does not predict the 
presence of organic hazes. Nor did previous photochemical models (Liang et al 2003, Liang 
et al 2004). Liang et al (2004) showed that simple hydrocarbons would not condense to 
form photochemical smogs in hot solar composition atmospheres. They concluded that 
photochemical smogs would not be present in the hottest hot Jupiter atmospheres that 
they addressed, but they did not consider cooler atmospheres where hydrocarbons would be 
more stable. 

Zahnle et al (2009) developed a new ID atmospheric chemistry model that addresses 
photochemical and thermochemical disequilibrium in hot Jupiter atmospheres. The focus 
of our first paper was on sulfur chemistry at temperatures between 1200 and 2000 K. Sulfur 
is a relatively abundant volatile element that equilibrium calculations suggest should be 
present as H 2 S (Visscher et al 2006). Our disequilibrium calculations suggest that sulfur 
should be present at observable altitudes as S 2 or HS. Sulfur is important to the chemistry 
because H 2 S is the most reactive primary molecule in solar composition atmospheres at 
temperatures between 800 K and 2000 K. Its thermal decomposition is the biggest source 
of free radicals and H in particular. We found that HS and S 2 should be abundant at 
the millibar level, and we suggested that most incident light from 250 nm to 450 nm 
would be absorbed by stratospheric S 2 and HS. The resultant stratospheric heating can be 
considerable in planets that are not hot enough for TiO and VO to evaporate in significant 
quantities. 

In passing we also confirmed an earlier suggestion (Fegley and Lodders 2002, Visscher 
et al 2006) that C0 2 should be a very good indicator of metallicity. We showed that this 
result applies when kinetics are included and that the C0 2 abundance and the C0 2 /CO 
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ratio are insensitive to T in the range 1200 < T < 2000 K. 



3. Model 

We use a standard ID kinetics code to simulate atmospheric chemistry of hot irradiated 
giant planets. Volume mixing ratios /$ of species i are obtained by solving continuity 



and force (flux) 



& = hJi ( ^ - ^) - fa + K ZZ N) ^ (2) 



equations. In these equations N is the total number density (cm -3 ); Pj — LiNfi represent 
chemical production and loss terms, respectively; <pi is the upward flux; bi a , the binary 
diffusion coefficient between i and the background atmosphere a, describes true molecular 
diffusion; m a and rrii are the molecular masses of a and i; and K zz , the eddy diffusion 
coefficient, parameterizes vertical mixing as diffusion. We have implemented molecular 
diffusion through H 2 by setting bi a = 6 x 10 19 (T/1400) ' 75 cm _1 s _1 — appropriate for CO 
- for all the heavy species. This is a reasonable choice for present purposes. The physical 
meaning of b ia is the ratio of the relative thermal velocities of two species to their collision 
cross section. In the present circumstances, the relative thermal velocity is effectively that 
of H 2 , so the only important source of variation in 6 ia is in the different diameters of the 
molecules. 

In practice Eqns [1] and [2] are solved as a system of second order partial differential 
equations for fi(z,t). Steady state solutions are found by integrating the equations through 
time (typically 10 9 years) using a fully implicit backward-difference method. Some aspects 
of the code are briefly described in Zahnle et al (2009). Other aspects are more completely 
described as applied to ancient Earth (Zahnle et al 2007) and Mars (Zahnle et al 2008). 
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Here we fully describe the chemical system. 

In this study we solve 528 chemical reactions and 33 photolysis reactions for 58 
chemical species. The 58 species are listed in Table 1. The reactions and references are 
listed, and where appropriate discussed, in Tables 2 and 3. In our scheme the hydrocarbon 
chemistry is truncated at C 2 H m (with the exceptions of C 4 H and C 4 H 2 ). This means that 
polymerization beyond C 2 H m is not included. Thus when conditions favor polymerization 
carbon pools in C 2 H m , because longer carbon chains are not allowed. 

New to this model are several hydrocarbon species added since Zahnle et al (2009). 
Methanol (CH 3 OH) and the radicals CH 3 and H 2 COH allow us to fully describe the 
known hydrogenation pathways from CO to CH 4 . These pathways need to be included 
because equilibrium chemistry predicts that CO should convert to CH 4 under warm 
Jupiter conditions. In practice these reactions are quite unimportant and CO is effectively 
indestructible under any conditions encountered in this study. Given that CO does not 
hydrogenate by gas phase reactions to any significant extent, we see no call at this time to 
include speculative pathways through the extremely unstable N 2 H radical that might allow 
hydrogenation of N 2 to NH 3 . 

In the original preprint edition of this paper we did not include addition reactions of 
hydrocarbons with OH, such as C 2 H 4 + OH + M — >■ adduct + M. The precise nature of the 
adduct is generally unknown, nor is it known in general how the adduct reacts. But as 
discussed above, the C-0 bond created cannot be easily undone in the atmosphere. In the 
preprint we considered only the competing H abstraction path C 2 H 4 + OH — >■ C 2 H 3 + H 2 0. 
This is a reasonable approximation at high temperatures and low pressures where the 
abstraction path is favored over the (3-body) addition. But by neglecting the reaction path 
in which the OH directly attacks a carbon atom, we underestimated the importance of 
oxidation at higher pressures and lower temperatures, and thus the preprint overstated the 
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stability of organic molecules. 

We have since upgraded the chemical model to explicitly include the adduct species 
C 2 H 2 OH, C 2 H 4 OH, and HCNOH. Reaction rates for making the adducts are described in 
the literature (R510, R512, R514). What happens to the adducts subsequently is neither 
simple nor fully known. For simplicity we have limited the options to reactions with atomic 
hydrogen, which is by far the most abundant free radical in these simulations, and we 
have chosen simple products that do not artificially reduce the number of free radicals 
(R511, R513, R515). The specific products for R511 and R515 were picked by analogy to 
the documented reaction of C 2 H 4 OH (R513). These improvements in the reaction scheme 
lead to moderate changes from the original preprint in the direction expected: in the new 
calculations the temperature threshold for hydrocarbon polymerization is about 100 K 
cooler. 

An important simplification in this study, as in its predecessors, is that we assume 
isothermal atmospheres with constant eddy diffusion coefficients. This choice facilitates 
presenting the results of parameter studies, as T and K zz prove to be the most important 
parameters. Published temperature-pressure profiles for hot jupiters are themselves models. 
These models depend explicitly on the assumed equilibrium chemical composition of the 
atmosphere; the absence (or presence) of arbitrary amounts of hazes or clouds; time of 
day or zonal circulation; etc. This study emphasizes the chemistry, and its intent is to 
provide the next step in the iteration between chemical and radiative-convective models. A 
"realistic" temperature profile requires at least two parameters (an effective temperature 
and a mean infrared opacity at minimum), and at least three more (an internal luminosity, 
visible opacity, and an insolation geometry) for full fidelity. Application to more realistic 
temperature profiles will be deferred to future work. 

What K zz (z) should be is not well constrained. Values ranging from 10 3 at the top of 
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the troposphere to 10 7 cm 2 /s at the top of the stratosphere seem to be useful for Jupiter. 
Because the amplitudes of upward propagating waves grow as p~ - 5 , it is often assumed 
that K zz oc p~ - 5 . The strong solar forcing of hot Jupiters can lead to much more energetic 
motions than seen on Jupiter itself. Showman et al (2009) developed 3D models of the 
circulation in hot Jupiters HD 189733b and HD 209458b. They suggested that "eddy 
diffusivities at 1 mbar are ~ 10 11 cm 2 sec 1 ." For simplicity we assume K zz . 

The background atmosphere is assumed 84% H 2 and 16% He. The relative abundances 
of C, N, O, and S are solar. The elements other than H and He are scaled as a group 
according to metallicity. Metallicity is varied from —0.7 < [M/H] < 1.7. Surface gravity is 
g — 20 m/s 2 . Surface gravity is important to optical depth and therefore important to what 
a planet looks like, but it has little effect on the chemistry (Zahnle et al 2009). Incident 
UV radiation is 100 x greater than at Earth and the solar zenith angle is 6 = 30° unless 
otherwise noted. When examining the effect of photolysis we vary the UV flux from 1 to 
1000 times solar. 

The upper boundary is a zero flux lid at ~ 1/zbar. For the lower boundary we assume 
thermochemical equilibrium mixing ratios of H, H 2 0, CO, CH 4 , NH 3 , N 2 , and H 2 S (Lodders 
and Fegley 2002, Visscher et al 2006). Other species are assumed to vanish at the lower 
boundary. This ensures that photochemical products flow into the deep atmosphere, where 
they are presumptively recycled. 

The lower boundary pressure was immaterial to our previous study because 
temperatures were very high and the relevant chemical reactions very fast. But at 
temperatures below 1000 K, the carbon and nitrogen chemistries become sensitive to the 
choice of lower boundary conditions. Here we impose chemical equilibrium at the lower 
boundary. The puzzle is where to put it. Chemical equilibrium is not actually expected in 
our model. Vertical mixing at 1000 K lifts stable molecules like NH 3 faster than it reacts 
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and, even at 100 bars, gas phase chemistry provides no effective means of converting CO to 
CH 4 and N 2 to NH 3 , and is only beginning to convert excess H to H 2 . In practice it is likely 
that CO is converted to CH 4 and N 2 is converted to NH 3 on grains, in particular Fe-Ni 
grains if present. Here we consider arbitrary lower boundary pressures of 1 and 100 bars. 

4. Results 

The lower temperatures considered in this study result in a much richer chemistry 
than we encountered in Zahnle et al (2009). Temperature, vertical mixing, UV and EUV 
irradiation, metallicity, and the choice of lower boundary conditions can all be important 
parameters in the 800-1200 K temperature range. We begin with an investigation of how 
our results vary with a single parameter when the other parameters are held fixed. Figures 
[TH3] are classic spaghetti plots used to show how computed atmospheric compositions vary 
with altitude. Figure [T] addresses temperature, Figure [2] addresses vertical mixing, Figure 
|3] addresses UV irradiation, and Figure ES addresses metallicity. Figure [1] is also used to 
illustrate the sensitivity to the lower boundary condition in one particular case, and Figure [2] 
is also used to compare computed disequilibrium chemical compositions to thermochemical 
equilibrium compositions at 1000 K. 

Spaghetti plots have their place, but it can be hard to see the entree for the noodles. 
To describe the bigger picture, we use pie charts that show computed chemical compositions 
at a certain key altitudes across the general parameter survey. The pies are arrayed in 
rows of constant K zz , according to temperature or metallicity. We use these to discuss the 
propensity of some atmospheres to form soots. The discussion begins with and concentrates 
on hydrocarbons, but we also provide brief overviews of nitrogen and sulfur chemistry. 
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4.1. Parameters 

Figure [1] illustrates vertical profiles of representative isothermal atmospheres at 1200, 
1000, and 800 K. These use K zz = 10 7 cm 2 /s for vertical mixing and [M/H] = 0.7 for 
metallicity, both comparable to those seen in Jupiter. The overall structure of these 
atmospheres resembles a poorly functioning gas stove, with methane introduced at the 
bottom and flowing upward, where it is either oxidized to CO or made into bigger organic 
molecules. At 1200 K, methane is quickly and cleanly oxidized, but at 800 K methane 
combustion generates mostly ethylene (C2H4) and HCN, to the point where ethylene 
becomes more abundant than methane. This does not mean that C2H4 actually would 
be more abundant than CH 4 (our model truncates polymerization at C^H^), but it does 
mean that polymerization is a major methane reaction pathway, and it suggests that 
hydrocarbons bigger than C2H4 are likely to form. 

The fourth panel of Figure [T] illustrates the effect of changing the lower boundary 
pressure. Placing the lower boundary at 1 bar puts more methane and less ammonia at 1 
bar than there would be if the boundary were at 100 bars. Ammonia is more abundant in 
the 100 bar simulations because (i) its equilibrium abundance is higher at higher pressure 
and (ii) its chemical destruction is kinetically inhibited at 1200 K and below. Hence the 
deeper lower boundary puts more NH3 in the stratosphere. Similar behavior would be seen 
in CH4 at T < 800 K. But at 1000 K, putting the lower boundary at 1 bar introduces an 
artificial source of methane that results in a big upward flux that in turn results in more 
ethylene. As both the 1 bar and 100 bar lower boundaries are arbitrary in this study it is 
not useful to take this discussion any further here. 

The chief factors behind the temperature dependence seen in Figure [1] are the high 
abundances of H and H2 and the different thermal stabilities of methane and water. Atomic 
H is more abundant than it would be in equilibrium. Disequilibrium of H and H2 is closely 
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tied to the disequilibria of H 2 S, HS, S, and S 2 . An important reaction that does not occur 
at significant rates is the 3-body recombination of H to H 2 : 

H + H + MoH 2 + M. Rl 

The computed H/H 2 ratio seen in the kinetics models is much higher than would be 
predicted by equilibrium. High disequilibrium H abundances were also a prominent feature 
in the models of Liang et al (2003) and Zahnle et al (2009). 

Both H 2 and CH 4 are attacked by H: 

H + CH 4 ->■ H 2 + CH 3 R45 

H + H 2 -> H 2 + OH. R6 

R6 is significantly endothermic, while R45 is nearly neutral, but both R45 and R6 have 
significant activation energies. Because R6 is endothermic, the temperature barrier is 
higher, so that water becomes unreactive at a higher temperature than methane, while the 
reverse of R6, 

H 2 + OH -> H + H 2 0. R5 

remains fast at 1000 K. The net effect of these trends is that, at 1200 K and above, water is 
very reactive and the gas is generally oxidizing, while methane reacts so quickly that it fails 
to reach altitudes significantly above the model's lower boundary. When the temperature is 
below 1000 K, water becomes stable enough that OH is scarce, while CH 4 becomes stable 
enough to be abundant in the atmosphere, yet unstable enough to be very reactive. As 
a consequence the effective C/O ratio in the reacting gases is high and the products are 
hydrocarbons. 

Mismatches between the kinetics model and thermochemical equilibrium at the lower 
boundary cause flows into or out of the lower atmosphere. The most important mismatches 
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stem from kinetic inhibition of the 3-body reactions that hydrogenate CO to CH4 and N2 
to NH 3 , and the 3-body reactions that reconstitute H 2 from H. The latter results in a high 
disequilibrium abundance of H. Very high abundances of H are a prominent feature of all 
the atmospheres in Figures [1] and [21 The overabundance of H pushes R45 farther to the 
right than it would be in thermochemical equilibrium. 

Figure [2] illustrates the influence of different rates of vertical transport in 1000 K 
atmospheres. The overall trend is that faster vertical transport lifts CH4 to higher altitudes, 
and thus raises the altitudes where C2H2, C2H4, and HCN form. Faster transport increases 
the efficiency of hydrocarbon formation vis a vis CO formation, as also seen in the CO 
profiles, and it changes the mix of products to favor C2H2 and HCN over C2H4. Apparently 
soot formation is likeliest if K zz > 10 10 cm 2 /s. 

The lower right-hand panel of Figure [2] shows what a 1000 K chemical equilibrium 
atmosphere would look like. The chemical equilibrium atmosphere can be thought of as 
having no vertical mixing (K zz = 0) and no photochemistry (1 = 0), although it would not 
be possible to reach equilibrium by setting K zz = in the kinetics model, because there are 
no efficient pathways for making NH 3 from N 2 or CH 4 from CO (in practice we found it 
difficult to obtain converged solutions with K zz < 10 5 cm 2 /s). In general, the disequilibrium 
kinetic atmospheres contain a much richer mix of chemical species. Product species are 
generally more abundant in the kinetics models, but so too is ammonia. On the other hand, 
at the altitudes where they would be observed, CO and C0 2 tend to be present at nearly 
their equilibrium abundances. 

Although the disequilibrium chemistry seen in our models is primarily due to thermal 
chemistry and vertical transport, photolysis can be important, especially at high altitudes 
and lower temperatures. Figure |3] compares the effects of photolysis at higher and lower 
levels of UV and EUV irradiation than used in Figures [1] and |2j Figure [3] shows that UV 
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and EUV photolysis inhibits hydrocarbon formation. Soot's precursors grow better in 
the dark. The consequences of photolysis in warm Jupiter atmospheres are more or less 
opposite to what is seen on Titan. In cold atmospheres like Titan's, methane photolysis 
leads to production of more complicated hydrocarbons. But in warm Jupiters, methane is 
broken up and hydrocarbons are made by thermochemistry. Photolysis's role is to free OH 
radicals from water. Oxidation by OH inhibits or prevents hydrocarbon polymerization. 
All of our atmospheres become more oxidized, and presumably less hazy, at high altitudes 
due to photolysis of H 2 0. In effect, the topmost atmosphere is bleached by photolysis. The 
chemocline between the oxidized upper atmosphere and the reduced lower atmosphere is 
especially pronounced in the 800 K atmosphere in Figure [1] where the temperature is low 
enough that thermal decomposition of water is relatively minor. 

Figure H] addresses metallicity. Enhanced metallicity favors molecules such as N 2 and 
CS2 that contain multiple elements other than H over molecules such as CH 4 and NH 3 that 
have only one. Enhanced metallicity also creates a more oxidized gas, because we have 
assumed a solar C/O ratio less than unity. These two trends work together to favor CO 
and especially CO2 at higher metallicities, but the two trends work to cross purposes with 
respect to hydrocarbon synthesis and soot formation, because the more oxidized conditions 
suppress hydrocarbon formation. 



4.2. Hydrocarbons 

Methane polymerization begins with R45, H + CH4 — > H 2 + CH3. Detailed analysis 
reveals that the chief pathway for making C 2 H n hydrocarbons is the reaction between CH3 
radicals to make ethane, 

CH 3 + CH 3 + M -> C 2 H 6 M. R289 
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This came as something of a surprise to us because ethane is not itself very abundant in 
our models. Evidently the more abundant C2H n hydrocarbons are formed by reactions such 
as C2H n + H — > H 2 + C2H n _ 1 , which are also pushed to the right by the high disequilibrium 
abundance of H. 

The influence of K zz on carbon chemistry is viewed from a different angle in Figure 
Here we show how carbon is allocated at the altitude that is most favorable to hydrocarbon 
production. To do this we divide carbon into 3 categories. The primary category is CH4, 
the source molecule. The other two are oxidized species containing CO bonds (CO, CO2) 
and the more reduced products that do not contain CO bonds (C2H 2 , C2H4, HCN, CS2, 
etc). Call these xcn 4 , %co, an d xq 2 , respectively. Total carbon is S c = £ch 4 + ^co + x c 2 - 
We define the most favorable altitude pf for organic synthesis to be where xc 2 /^c reaches 
its maximum value. Whilst every C-0 bond is doomed to end in the irreversible production 
of CO, it is possible that much of the C in xq 2 can end in larger molecules such as PAHs or 
in the irreversible production of refractory soots. 

Figure M generalizes these results to different atmospheres as pie charts on T-K zz and 
[M/H]-iT zz grids. The relative amounts of carbon in different species is presented at the 
altitude pj. It is an optimistic figure, in the sense that it shows the most favorable conditions 
for carbon polymerization in each model. Ethylene is favored by cooler atmospheres with 
strong vertical mixing, while CO and CO2 are favored by higher temperatures, higher 
metallicity, and weaker vertical mixing. Acetylene — soot's precursor — is favored by very 
strong vertical mixing. The dependence of C2H2 itself on temperature is relatively weak, 
but the dependence of CO on temperature suggests that competition with oxidation makes 
higher temperatures less suited to soot formation. Figure [6] also shows a clear negative 
dependence on metallicity that suggests that low metallicities are more favorable to soot 
formation. The origin of the effect is that metals are on balance oxidizing because the C/O 
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ratio in a solar gas is less than one. Therefore adding more metals creates more oxidized 
conditions, other things equal. 

Ethylene, acetylene, and hydrogen cyanide are first generation products of methane. 
Benzene (CeH 6 ) would be at best a third generation product (e.g., three acetylenes), and 
PAHs like coronene (C24H 12 ) and ovalene (C32H44) are probably two or three generations 
further evolved. At each step some carbon is lost to CO. If it takes n generations to make 
something resembling a soot, we might expect the efficiency of soot formation to go like 
Xq 2 . This is greater than 10% for xc 2 > 0.8 and n = 10, so we infer that atmospheres where 
Xc 2 > 0.8 have potential to develop significant organic hazes. 

PAHs interact strongly with visible light (Lou Allamandola, pers. comm.). Coronene is 
golden yellow, ovalene is red, and still bigger ones tend to black. Thus PAHs by themselves 
may suffice to explain high altitude optical absorption in hot Jupiters; it may not be 
necessary that they agglomerate into grains, albeit grains are likely to be more refractory. 

Figure [7| uses pie charts to provide an overview of the hydrocarbon chemistry that 
is much easier to interpret than the spaghetti plots shown in Figures (THU (which are 
informative but daunting). Here the pie charts show carbon speciation at two key altitudes, 
chosen to approximate the levels where the optical depths would be of order 0.01 and 1 
in clear air. Optical depths r = 0.01 and r = 1 suggest what might be seen in grazing 
incidence (primary eclipse) and in direct imaging (secondary eclipse), respectively. To 
estimate the optical depths we have approximated the pressure-dependent Rosseland 
infrared mean opacities from Freedman et al (200x), for equilibrium chemistry for clear 
air of solar composition, by k ~ 0.03p 0,6 cm _2 /g, where p is in bars. This is a very crude 
fit, and the assumptions on which it is founded are of questionable relevance, yet the 
results may have some value for illustrative purposes because the strongest infrared opacity 
source would in all cases be water vapor. Pressure levels corresponding to r ~ 1 and 
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r m 0.01 for four different metallicities are listed in Table 4. Perhaps the most obvious and 
robust conclusion to draw from Figure [7J (and the analogous figures for nitrogen and sulfur 
discussed below) is that the apparent chemical composition of a planet depends on how it is 
viewed (Fortney 2005). In Figure [7J methane is much less prominent in transit observations 
than it is in observations of the photosphere. 

4.3. Nitrogen 

Figure |8] provides an overview of nitrogen speciation that is analogous to Figure [7J 
for carbon. Ammonia is favored by low temperatures, low metallicity and strong vertical 
mixing, and N2 is favored by the opposites. HCN, which forms most efficiently when 
ammonia is present and methane is polymerizing, is most abundant on the boundary 
between NH 3 and N 2 . The chief pathway for making HCN in our models is 

NH 2 + C 2 H 3 -> CH 3 + HCN. R337 

This reaction is known to make an adduct; we have assumed that the C-N bond in the 
adduct leads inevitably to HCN, much as we have assumed that the C-0 bond in adducts 
made by OH and C2H4 or C2H 6 lead inevitably to CO. That a poorly characterized reaction 
like R337 should play a key role is not ideal. Analogy to R289 suggests that the true path 
to HCN should pass through methylamine, CH3 + NH 2 — > CH3NH2. We have not yet 
included this pathway in our chemical scheme because we had no inkling that it should be 
important. This means that our predictions for HCN abundances are more uncertain than 
those for the hydrocarbons, and that our HCN abundances are likelier to be underestimated 
than overestimated. 
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4.4. Sulfur 

Sulfur species are useful because they interact strongly with blue, violet, and UV 
light, which makes them potentially observable, especially in transit. Sulfur speciation 
is illustrated in Figures and [TQJ The former uses spaghetti plots to show the vertical 
structure of the sulfur species in a few selected models. The latter is an overview chart 
analogous to Figures [7] for carbon and [8] for nitrogen. Equilibrium calculations predict 
that H 2 S would be the most abundant S-containing species under nearly any condition 
encountered in these atmospheres, as is seen above in Figure |2j 

Zahnle et al (2009) discussed sulfur chemistry at higher temperatures (T > 1200 
K). At these temperatures S 2 and HS are important at high altitudes but neither CS 
nor CS2 were important. In the cooler atmospheres investigated here both CS and CS2 
become important in parallel to hydrocarbons like ethylene, and for the same reasons: 
the increased availability of reactive methane, and the decreased availability of OH from 
water. Consequently cooler atmospheres with strong vertical mixing favor disequilibrium 
production of CS and CS2. Diatomic sulfur (S2) is predicted to be most abundant at high 
altitudes at higher temperatures, higher metallicities, and weaker vertical mixing. These 
dependencies are obvious in Figure [101 It is interesting that S 2 and acetylene (soot) appear 
to be mutually exclusive, the former indicating relatively metal-rich placid conditions, the 
latter indicating strong vertical mixing and low metallicity. 

5. Discussion 

The chemistry that we have modeled is unremarkable. Most of the important reactions 
involving small organic molecules and radicals have been studied both experimentally and 
theoretically. Thus the prediction that methane can react to make ethylene rather than CO 
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seems secure. The presumption that this will lead to bigger hydrocarbons than ethylene 
seems equally secure. Whether this process continues up to PAHs and soot is less secure. 

On the surface our results appear to disagree with what Liang et al (2004) found 
previously. However, there are enough differences between their simulations and ours that 
the apparent disagreement is misleading. The higher hydrocarbons that form in our models 
are primarily products of disequilibrium thermochemistry, not photochemistry. They form 
because recombination of H to H 2 is kinetically inhibited and consequently the atomic 
hydrogen density is much higher than it would be in equilibrium. Liang et al find less H 
than we do because they did not include sulfur. The most important reactions governing 
the H abundance are those involving sulfur, especially 

H + H 2 S H 2 + HS R156 

H 2 + HS H 2 S + H R157 

and 

S + H 2 <-> H + HS R168 

H + HS <-> H 2 + S. R169 

Reactions R156 and R157 are generally the most frequently occurring reactions in the 
atmosphere by a wide margin (Zahnle et al 2009). Reactions R168 and R169 are most 
important at higher altitudes (lower pressures). Photolysis has little to do with these. 

Another difference between our study and Liang et al (2003, 2004) is that Liang et al 
attempt to use realistic approximations to T(z) and K zz (z). Constant temperature and 
constant K zz is a feature of our study. Our emphasis here is on the survey of parameter 
space. Temperature and vertical mixing are the two key parameters. The cleanest way 
to look at the problem is therefore to use T and K zz as the independent variables. More 
complicated atmospheric profiles of either T(z) or K zz {z) requires introducing additional 
parameters, which introduces an additional level of complexity that may not be justified. 
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For K zz , Liang et al use K zz oc N~ os . Their assumption gives K zz = 2.4 x 10 7 cm 2 /s at 
1 bar, 1.5 x 10 9 cm 2 /s at 1 mbar, and 1.0 x 10 11 cm 2 /s at 1 /ibar (all at 1200 K). This range 
is similar to our high K zz cases. For T(z), Liang et al use three specific profiles derived 
from detailed radiative transfer model results. All three are very hot at 1 bar (the coolest 
is 1800 K). These temperatures are much higher than anything we address here. In our 
model methane is immediately destroyed if T > 1200 K. We agree with Liang et al (2004) 
that very hot planets ought not to have hydrocarbon hazes. 

As a test we computed chemistry using the nominal p-T profile from Fortney et al's 
(2009) most recent model of HD 189733b. Their model is hot (T « 1500 K) below 3 bars 
and cool (T m 900 K) above 10 mbars, with a smooth transition in between. We note that 
Fortney et al include only Rayleigh scattering, so that their computed albedo is relatively 
low, and thus their T(z) may not closely resemble that of a hazy planet. For the simulation 
we set K zz = 10 8 . We find very little methane above the lower boundary and no significant 
polymerization. We then reduced the lower temperature to 1200 K, and then to 1100 K. 
At 1200 K, NH3 is abundant but hydrocarbons are not, while at 1100 K we see rampant 
hydrocarbon polymerization. These results are quite consistent with Figure |6j 

In passing, we note that solar C/O is also an assumption. Snowline models of giant 
planet formation can give lower C/O, while the less well known tar-line model (Lodders 
2004) predicts higher C/O. The true C/O ratio of Jupiter is not yet known. High C/O 
changes the chemistry dramatically; among other things, soots would be strongly favored at 
higher temperatures. The possibility that a planet like HD 189733b might have high C/O 
should be in play. 
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6. Conclusions 

We find that vertical mixing in hot J up iters with temperatures below 1000 K can 
provide environments suitable for turning methane into ethylene, acetylene, and hydrogen 
cyanide. We speculate that such atmospheres would be conducive to generating PAHs and 
perhaps soots. We also reiterate the point that planetary compositions inferred from transit 
observations, which probe high altitudes (low pressures), can differ markedly from those 
inferred from reflected or emitted light from the same planet. In general, thermal chemistry 
and vertical mixing are more important to disequilibrium chemistry in warm or hot Jupiters 
than is photochemistry. Photochemistry plays a bigger role at low temperatures. 

The presence of abundant PAHs or soots at high altitudes could account for the 
up-to-now unique properties of HD 189733b, if HD 189733b proves to be cool enough, or 
its metallicity low enough, or its C/O ratio high enough, or its vertical mixing as vigorous 
as Showman et al (2009) suggest. This would make HD 189733b the first of a new class of 
warm Jupiters. As noted, Swain et al (2009) reported that HD 189733b's metallicity might 
be rather low. On the other hand HD 189733 is a relatively strong source of UV (Knutson 
2010), which makes photolysis rates high, a factor that we have shown will suppress high 
altitude soot formation. Whether our thermometric interpretation of HD 189733b actually 
applies to HD 189733b (recent models suggest that the planet is too hot if it has solar C/O) 
is not the main point: many organic- rich hazy warm Jupiters are sure to be discovered in 
the near future. 
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Table 1 

Chemical Species 

Elements Species 

hTo h, h 2 , o, oh, H 2 0, 2 , 0( 1 D) 

CO CO, C0 2 , HCO, H 2 CO, H 2 COH, CH 3 0, CH 3 OH 

CH n C, CH, CH 2 , 1 CH 2 , CH3, CH4 

C m H„ C 2 , C 2 H, C 2 H 2 , C 2 H 3 , C 2 H 4 , C 2 H 5 , C 2 H 6 , C 2 H 2 OH, C 2 H 4 OH, C4H, C4H 2 

N N, N 2 , NH, NH 2 , NH 3 , CN, HCN, H 2 CN, HCNOH, NO, NS 

S S, HS, H 2 S, S 2 , S 3 , S4, Sg, Sg b , CS, CS 2 , OCS, HCS, H 2 CS, SO, HSO, S0 2 

a - Ring (ground state) 

b - Linear (excited state) 
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Table 2 

Chemical Reactions 





Reactants a 




Products 


Rate" 






Reference 


Rl 


H + H + M 


-4 


H 2 4- M 


8.8x10" 


" 33 (T/298)" u - bU 




Ba92 




H + H 


— >■ 


H 2 + M 


2.0x10" 


-10 






R2 


H 2 + M 


-4 


H + H +M 


1.5x10" 


-09 e -48400/T 




Ba92 


R3 


H 2 + 


-4 


OH + H 


3.5x10" 


13 (T/298) 2 ' 67 e" 


-3160/T 


Ba92 


R4 


H + OH 


-4 


H 2 + 


1.7x10" 


13 (T/298) 2 ' 64 e" 


-2240/T 


rev3 


R5 


H 2 + OH 


-4 


H 2 + H 


1.6x10" 


12 (T/298) L6 V 


-1660/T 


Ba92 


R6 


H9O + H 


— y 


H9 + OH 


6.9 x 10" 


12 (T/298) 1 - 60 e 


-9270/T 


Ba92 


R7 

-L L 1 


OH 4- OH 


— V 

7 


HoO 4- 


1 7x10" 


~ 12 (T/298) 1 - 14 e 


-50/T 


Ba92 


R8 


H + 


— y 


OH 4- OH 


1.8 x 10" 


11 (T/298) ' 95 e 


-8570/T 


Li91 


R9 


O + H + M 


— y 


OH 4- M 


4.3 x 10" 


" 32 (T/298)~ im 




Ts86 




+ H 


— y 


OH 4- M 


2.0 x 10" 


-11 






RIO 


OH + M 


— y 


+ H 4-M 


4.0 x 10" 


-09g-50000/T 




Ts86 


Rll 


H 4- 


— V 

7 


OH 4- hv 


9 9x10" 


-19 fT/298r ' 38 




Mi97 


R12 


4- 4- M 


— \ 

7 


4 M 


1.0 x 10" 


" 33 fT/298r L0 ° 




Wa84 

V V OjUt: 




4- 

\y vy 


— \ 

7 


Oo 4- M 


2. Ox 10" 


-11 






R13 


Oo 4- M 


— y 


4- O 4-M 


1.0 x 10" 


" 08 (T/298)' im 


e -59400/T 


Ts86 


R14 


+ OH 


— y 


O9 + H 


2.4 x 10" 


-ll e -353/T 




Ba92 


R15 


H + 2 


-4 


OH + 


3.3x10" 


-10 e -8460/T 




Ba92 


R16 


H 2 + 2 


-4 


OH + OH 


4.0x10" 


- n (T/298) a4 V 


-35100/T 


Ka05 


R17 


CO + M 


-4 


C 4- 4-M 


2.7x10" 


"° 3 (T/298r 3 ' 52 


gl30000/T 


Ba92 


R18 


c + 


-4 


CO + M 













C + 


-4 




n 








R19 


C + OH 


-4 


CO + H 


1.1x10" 


" 10 (T/298) - 50 




Mi97 


R20 


C + OH 


-4 


CH4O 


2.2x10" 


- 11 (T/298) - 50 e" 


-14800/T 


Mi97 


R21 


C + 2 


-4 


CO + 


1.6x10" 


-11 




Ba92 


R22 


H + CO 


-4 


C + OH 


2.8x10" 


09 (T/298) a50 e 


-76900/T 


rev21 


R23 


CO + OH 


-4 


C0 2 4- H 


1.8x10" 


" 14 (T/298) L89 e 


583/T 


Li07 


R24 


C0 2 + H 


-4 


CO + OH 


6.1x10" 


11 (T/298) - 64 e 


-12500/T 


rev23 


R527 


C 2 H 5 + H 2 CN -4 


HCN + C 2 H 6 


7.7x10" 


-12 






R528 


NH 2 + C 2 H 5 


-4 


NH + C 2 H 6 


3.0x10" 


" 13 (T/298) 1 - 00 e 


-4400/T 


Xu99 



a — M refers to the background atmosphere, principally H 2 and He; units of density [cm 3 ]. 



b — 2-body reaction rates are in cm 3 s 1 ; 3-body rates are in cm 6 s 
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Table 3 

Photolysis Reactions 





Species 




Products 


Rate a 


Reference 


PI 


H 2 + hz/ 


->■ 


OH + H 


1.6 xlCT 03 


Sa03 


P2 


C0 2 + his 


->■ 


CO + 


2.7xl0-° 7 


Ok78,Hu92 


P3 


C0 2 + hz/ 


->■ 


CO + 0( 1 D) 


3.0 xlO" 05 


Ok78,Hu92 


P4 


2 + hi/ 


->■ 


+ 


9.6 xKT 06 


Sa03 


P5 


2 + hi/ 


->■ 


+ 0( 1 D) 


5.1 xlO^ 04 


Sa03 


P6 


NO + hi/ 


->■ 


N + 


3.7xl0-° 6 


Sa03 


P7 


H 2 S + hi/ 


->■ 


HS + H 


2.7 xl(T 02 


Ok78,Hu92 


P8 


NH 3 + hi/ 


->■ 


NH 2 + H 


1.1 xlO -02 


Ok78,Hu92 


P9 


NH 3 + hi/ 


->■ 


NH + H 2 


5.5 xlO" 03 


Ok78,Hu92 


P10 


CH 4 + hi/ 


->■ 


CH 2 + H 2 


3.9 xlO^ 04 


Hu92 


Pll 


CH 4 + hi/ 


->■ 


CH 3 + H 


3.9 xlO" 04 


Hu92 


P12 


S0 2 + hi/ 


->■ 


SO + 


1.7xl0 -02 


Ok78,Hu92 


P13 


S0 2 + hi/ 


->■ 


s + 2 


6.9 xlO" 04 


Ok78,Hu92 


P14 


S0 2 + hi/ 


->■ 


S + +0 


1.7xl0 -05 


Ok78,Hu92 


P15 


SO + hi/ 


->■ 


S + 


4.4 xlO^ 02 


Ok78,Hu92 


P16 


CS 2 + hi/ 


->■ 


cs + s 


3.9 xlO" 01 


Mo81,Ah92 


P17 


OCS + hi/ 


— > 


CO + s 


2.5 xlO -03 


Sa03 


P18 


S 2 + hi/ 


->■ 


s + s 


1.7xl0 -01 


Za09 


P19 


S 3 + hi/ 


->■ 


s 2 + s 


l.lxl0+ 02 


Za09 


P20 


S 4 + hz/ 


->■ 


s 3 + s 


l.lxlO+ 01 


Za09 


P21 


S 8 + hz/ 


->■ 


s* 8 


1.7xl0 +0 ° 


Ka89 


P22 


S* + hz/ 


->■ 


S 4 + S 4 


3.5xl0 +0 ° 


Ka89 


P23 


C 2 H 2 + hz/ 


->■ 


C 2 H + H 


2.4 xlO" 03 


Ok78,Hu92 


P24 


C 2 H 4 + hz/ 


->■ 


C 2 H 3 + H 


5.3x10 UcS 


Ok78,H u 92 


P25 


C 2 H 6 + hz/ 


->■ 


C 2 H 5 + H 


1.1 xlO" 03 


Ok78,Hu92 


P26 


C 4 H 2 + hz/ 


— >■ 


C 4 H + H 


4.8 xlO" 03 


note 


P27 


H 2 CO + hz/ 




CO + H 2 


6.1 xl0~ 03 


Sa03 


P28 


H 2 CO + hz/ 


->■ 


HCO + H 


7.1 xlO -03 


Sa03 


P29 


H 2 CO + hz/ 


->■ 


CO + H +H 





Sa03 


P30 


HCN + hz/ 


->■ 


CN + H 


2.5 xlO" 03 


Hu92 


P31 


HSO + hz/ 


->■ 


HS + 


6.1 xl0~ 02 


note 


P32 


HS + hz/ 


->■ 


H + S 


1.0xlO+ 01 


Za09 


P33 


CH 4 + hz/ 


->■ 


CH + H 2 +H 


1.9 xlO" 04 


Ok78,Hu92 



a Photolysis rates at the top of the atmosphere for S uv = 100 and 30° zenith angle. 
P21. Photolysis of S 8 is presumed to create the linear radical S* 8 . This is a placeholder at high 
temperatures where S§ is not expected and at low temperatures where Sg condenses. 
P26. Assumes twice the cross section of C 2 H 2 . 
P31. Assumes cross section of H0 2 . 



Table 4 

Reference Pressures 



r = 0.01 r = l 



[M/H] p [mbars] p [mbars] 



0.2 


30 


500 


1 


10 


200 
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4 


70 


50 


1 
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Fig. 1. — Vertical profiles of selected H, C, N, O species in isothermal atmospheres at three 
temperatures (T = 1200, 1000, and 800 K). Fixed model parameters are K zz = 10 7 cm 2 /s, 
planetary ([M/H] = 0.7) metallicity, I = 100, and g = 20 m/s 2 . The lower boundary condition 
is thermochemical equilibrium at 100 bars. The fourth panel, to be compared to the panel 
above it, raises the thermochemical equilibrium lower boundary to 1 bar. At 1200 K methane 
is oxidized with little buildup of hydrocarbons. At 800 K methane's destruction leads to 
significant disequilibrium production of organic molecules, especially ethylene (C 2 H 4 ). The 
horizontal gray bars denote the most favorable altitudes for hydrocarbon formation pj as 
defined in the text. The decreasing abundance of stable species at high altitudes is caused 
by molecular diffusion. 
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Fig. 2. — Vertical profiles of selected species in isothermal 1000 K atmospheres and 
[M/H] = 0.7 as a function of vertical mixing. Strong vertical mixing — K zz > 10 9 — favors 
disequilibrium production of C2H4, HCN, and especially C2H2 (acetylene, the canonical soot 
precursor). The horizontal gray bars denote the most favorable altitudes for hydrocarbon 
formation pf as defined in the text. The fourth panel shows thermochemical equilibrium 
abundances at all altitudes at T = 1000 and [M/H] = 0.7. In general, the kinetics calcula- 
tions predict a much richer variety of species than equilibrium. 
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Fig. 3. — Vertical profiles of selected H, C, N, O species in isothermal atmospheres at 
higher and lower levels of UV and EUV irradiation. Fixed model parameters are T = 
1000 K, [M/H] = 0.7, and g = 20 m/s 2 . At 1000 K, photolysis is most important when 
vertical mixing is vigorous (K zz = 10 11 cm 2 /s). Photolysis is less important at K zz = 
10 9 , and at K zz = 10 7 (not shown) is important only for NH 3 and HCN. In general, UV 
is hostile to NH 3 and small hydrocarbons other than CH 4 , and through H 2 photolysis 
creates a more oxidizing environment. These effects suppress hydrocarbon formation and 
relegate hydrocarbon formation to a deeper level in the atmosphere. Soot's precursors, in 
particular, are strongly favored by darkness. The horizontal gray bars denote the most 
favorable altitudes for hydrocarbon formation pf. For 5=1 and K zz = 10 11 , this would lie 
above the top of the model. 
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Fig. 4. — Vertical profiles of selected H, C, N, O species in isothermal atmospheres at four 
metallicities. Relative abundances of C, N, S, and O are solar. Fixed model parameters 
are K zz = 10 9 cm 2 /s, T = 1000 K, S uv = 100, and g = 20 m/s 2 . The general trend is 
that higher metallicities result in more strongly oxidized mixtures. CO and especially CO2 
are good indicators of metallicity. The propensity to form soot appears to be similar for 
—0.7 < [M/H] < 0.7, as judged by C4H 2 abundances, but higher metallicity is unfavorable. 
The horizontal gray bars denote the most favorable altitudes for hydrocarbon formation pf. 
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Fig. 5. — Methane's fates at the most favorable altitude pf for soot formation as a function 
of vertical mixing K zz . Fates are shown at T = 1000 K (solid lines) and T = 900 K 
(dashed lines). Fixed model parameters are [M/H] = 0.7, S uv = 100, po = 100 bars, and 
g = 20 m/s 2 . Irreversible oxidization (xco) is denoted "CO and C0 2 ", interesting products 
(xc 2 ) are denoted "SC2H n ,HCN etc". High K zz favors organic molecules. At 1000 K, soot 
formation seems reasonably likely for K zz > 10 8 cm 2 /s. The right hand axis shows that pj is 
a monotonically increasing function of K zz . 
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Fig. 6. — Generalization of Figure [5] to other temperatures (left, [M/H] = 0.7) and metallic- 
ities (right, T = 1000 K). Common model parameters are S uv = 100, po = 100 bars, g = 20 
m/s 2 , and # = 30° incidence angle for UV radiation. Each pie-chart shows how carbon is 
allocated at pf for a specific model. To first approximation higher values of K zz map to 
higher altitudes (lower pf). Left Low temperatures and high K zz favor organic molecules. 
Very high values of K zz > 10 10 cm 2 /s mix CH4 to altitudes that are more favorable to C2H2, 
the traditional soot precursor. Right Low metallicity favors hydrocarbons generally and 
C2H2 and HCN particularly. High metallicity is quite oxidized. The apparent preference for 
CH 4 at high metallicity and low K zz is an artifact of there being essentially no hydrocarbon 
synthesis and pf being near the lower boundary. 
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Fig. 7. — Carbon speciation at r fa 0.01 (top) and r ~ 1 (bottom). These approximate what 
might be seen in grazing incidence in transit and in differential eclipses, respectively. For 
hydrocarbons vs. temperature, fixed parameters are [M/H] = 0.7, S uv = 100, and g = 20 
m/s 2 . For hydrocarbons vs. metallicity, fixed parameters are T = 1000 K, S uv = 100, and 
g = 20 m/s 2 . 
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Fig. 8. — Nitrogen speciation at r ^ 0.01 (top) and r « 1 (bottom). Cases and meanings 
are the same as for Figure [7J Ammonia is favored by low temperatures, low metallicity 
and strong vertical mixing, and N 2 is favored by the opposites. HCN, which forms most 
efficiently when ammonia is present and methane is polymerizing, is most abundant on the 
boundary between NH3 and N 2 . 
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Fig. 9. — Vertical profiles of sulfur species in isothermal atmospheres at different temper- 
atures and vertical mixing. S2 is favored by high temperatures and weak vertical mixing, 
while CS and CS2 are favored by lower temperatures and strong vertical mixing. H2S is 
the stable gas in the deep atmosphere. CS and CS2 are found in environments that also 
favor C2H4. S2 tends to form at the chemical boundary between the more reduced lower 
atmosphere and the more oxidized upper atmosphere. 
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Fig. 10. — Sulfur speciation at r « 0.01 (top) and r fa 1 (bottom). Cases and meanings are 
the same as for Figure [3 CS and CS2 are favored by cooler atmospheres with strong vertical 
mixing, while S2 is favored at high altitudes by higher temperatures, higher metallicities, 
and weaker vertical mixing. 



